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Improved Discrete Approximation of Laplacian of Gaussian 

This method reduces the amount of circuitry needed for filtering of video data. 

Lyndon B. Johnson Space Center, Houston, Texas 


An improved method of computing a 
discrete approximation of the Laplacian 
of a Gaussian convolution of an image 
has been devised. The primary advan- 
tage of the method is that without sub- 
stantially degrading the accuracy of the 
end result, it reduces the amount of in- 
formation that must be processed and 
thus reduces the amount of circuitry 
needed to perform the Laplacian-of- 
Gaussian (LOG) operation. 

Some background information is nec- 
essary to place the method in context. 
The method is intended for application 
to the LOG part of a process of real-time 
digital filtering of digitized video data 
that represent brightnesses in pixels in a 
square array. The particular filtering 
process of interest is one that converts 
pixel brightnesses to binary form, 
thereby reducing the amount of infor- 
mation that must be performed in sub- 
sequent correlation processing (e.g., 
correlations between images in a stereo- 
scopic pair for determining distances or 
correlations between successive frames 
of the same image for detecting mo- 
tions). The Laplacian is often included 
in the filtering process because it em- 
phasizes edges and textures, while the 
Gaussian is often included because it 
smooths out noise that might not be 
consistent between left and right images 


or between successive frames of the 
same image. 

For the purpose of processing digi- 
tized image data, the Gaussian and 
Laplacian values of a pixel of interest are 
approximated as weighted sums of 
brightnesses of the pixels in a square 
subarray (typically 3x3) of pixels cen- 
tered on the pixel of interest. The 
weights are represented by coefficient 
matrices, the elements of which corre- 
spond to pixels in the subarray. For ex- 
ample, 
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gives a reasonable approximation of the 
Laplacian for a 3 x 3 subarray. 

A typical prior state-of-the-art LOG algo- 
rithm operates in a sequential raster-ori- 
ented pixel stream, and the kernel is fac- 
tored as much as possible into 1x3 and 3 
x 1 components. To apply a 1 x 3 or a 3 x 
3 kernel to the pixel in a given row and 
column, it is necessary to retain the pixel 
stream in two raster-length delay lines so 
that the pixels in the adjacent rows and 
columns remain available for the compu- 
tation. Heretofore, 12 bits of precision 
have been needed to maintain accuracy 
sufficient for reasonable convolution re- 


sults through several stages. Usually, the 
Laplacian operation is performed first to 
normalize the image data about 0. This 
concludes the background information. 

In the present method, intermediate 
results are approximated in such a way 
that only 6 or possibly even as few as four 
bits are retained, yet the final result is 
still reasonable. If only 6 bits of precision 
are needed rather than 12, the size of 
the memory circuitry needed to imple- 
ment the delay lines can be halved. 
Thus, it should be possible to build 
smaller, lower-power filtering circuits. 

Heretofore, it has been common prac- 
tice in limited-precision arithmetic cir- 
cuitry to approximate large values by 
truncating them, eliminating the least 
significant bits. However, a detailed 
analysis of the arithmetic process shows 
that eliminating the bits of lowest order 
can lead to errors in the final conversion 
to binary representation. In the present 
method, the least significant bits are re- 
tained and large values are approxi- 
mated by a saturation technique based 
on the unconventional approach of dis- 
carding the highest-order bits. If the 
magnitude of a pixel value is larger than 
the largest magnitude that can be repre- 
sented in the result, then the pixel value 
is replaced by a value of the same sign 
(positive or negative) and the largest 
representable magnitude. 

The figure depicts an example of a cir- 
cuit that utilizes 2’s-complement encod- 
ing of negative numbers and that imple- 
ments saturation of a six-bit quantity to a 
four-bit quantity. The exclusive-OR gates 
examine the number of high-order bits 
to be eliminated, plus one, to determine 
whether they are already all alike. If the 
examined input bits are not all alike, a 
maximum positive or negative quantity is 
derived from the sign bit and gated 
through the multiplexer. This is easily ex- 
tended to handle a reduction of any 
number of bits. 

This work was done by Robert L. Shuler, Jr., 

of Johnson Space Center. Further infor- 
mation is contained in a TSP ( see page 1 ). 
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The Input Bits Are Passed through to the multiplexer as long as the magnitude of the input quantity 
does not exceed a preset four-bit saturation value; if it does exceed that value, then the output has 
the saturation magnitude and the same sign as that of the input. 
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